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Abstract 

We present a phase-space analysis of cosmology containing multiple scalar fields with posi- 
tive and negative exponential potentials. We show that there exist power-law multi-kinetic- 

K*- ■ potential scaling solutions for sufficiently fiat positive potentials or steep negative poten- 
00 ' 

Tjj- ■ tials. The former is the unique late-time attractor and the well-known assisted infiationary 

Q ' solution, but the later is never unstable in an expanding universe. Moreover, for steep neg- 

O ' ative potentials there exist a kinetic-dominated regime in which each solution is a late-time 

2 ■ attractor. We briefiy discuss the physical consequences of these results. 
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Scalar field cosmological models are of great importance in modern cosmology. The dark 
energy is attributed to the dynamics of a scalar field, which convincingly realizes the goal 
of explaining present-day cosmic acceleration generically using only attractor solutions . 
A scalar field can drive an accelerated expansion and thus provides possible models for 
cosmological inflation in the early universe [2]. In particular, there have been a number 
of studies of spatially homogeneous scalar field cosmological models with an exponential 
potential. There are already known to have interesting properties; for example, if one 
has a universe containing a perfect fluid and such a scalar field, then for a wide range of 
parameters the scalar field mimics the perfect fluid, adopting its equation of state jHI. These 
scaling solutions are attractors at late times |3] . The inflation [31 E] and other cosmological 
effect [Zj of multiple scalar fields have also been considered. 

The scale-invariant form makes the exponential potential particularly simple to study 
analytically. There are well-known exact solutions corresponding to power-law solutions 
for the cosmological scale factor a oc in a spatially flat Friedmann-Robertson- Walker 
(FRW) model [S", but more generally the coupled Einstein-Klein- Gordon equations for a 
single field can be reduced to a one-dimensional system which makes it particularly suited 
to a qualitative analysis [HlCIIj. Recently, adopting a system of dimensionless dynamical 
variables the cosmological scaling solutions with positive and negative exponentials 
has been studied ^2]- Usually there are many scalar fields with exponential potentials in 
supergravity, superstring and the generalized Einstein theories, thus multi potentials may 
be more important. In this paper, we will consider multiple scalar fields with positive and 
negative exponential potentials. We have assumed that there is no direct coupling between 
the exponential potentials. The only interaction is gravitational. A phase-space analysis of 
the spatially flat FRW models shows that there exist cosmological scaling solutions which 
are the unique late-time attractors, and successful inflationary solutions which are driven 
by multiple scalar fields with a wide range of each potential slope parameter A. 

We start with more general model with m scalar fields (pi, in which each has an identical- 
slope potential 



where = SttGn is the gravitational coupling and A is a dimensionless constant charac- 
terising the slope of the potential. Note that there is no direct coupling of the fields, which 
influence each other only via their effect on the expansion. The evolution equation of each 
scalar field for a spatially flat FRW model with Hubble parameter H is 
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subject to the Friedmann constraint 
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Defining 2m dimensionless variables 



x- = ^ y = ^ (4) 



the evolution equations can be written as an autonomous system: 
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x!- = -3x,\l-J2x^±X^-y^ (5) 



m 
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y'i = Vi \ - I 



where a prime denotes a derivative with respect to the logarithm of the scalar factor, 
= Ina, and the constraint equation Q becomes 



Throughout we will use upper/lower signs to denote the two distinct cases of ±Vi > 0. xj 
measures the contribution to the expansion due to the field's kinetic energy density, while 
±yf represents the contribution of the potential energy. We will restrict our discussion of 
the existence and stability of critical points to expanding universes with H > 0, i.e., y >0, 
and A > 0. Critical points correspond to fixed points where x[ = 0, = 0, and there are 
self-similar solutions with 

H 

-77-2= -^T.^" (8) 

i=i 

This corresponds to an expanding universe with a scale factor a{t) given by a oc or a 
contracting universe with a scalar factor given by a oc {—ty, where 

1 

P = oy^rn o (9) 

The system (jS} and © has at most one m-dimensional sphere embedded in 2m-dimensional 
phase-space corresponding to kinetic-dominated solutions, and (2™ — 1) fixed points, one 
of which is a m-kinetic-potential scaling solution. 

In order to analysis the stability of the critical points, we only consider the cosmologies 
containing two scalar fields. There are one unit circle 5* and three fixed points Ai, A2 and 
B listed in Table 1. Using the Friedmann constraint equation (jZj), we reduce Eqs.dS)) and 
(P) to three independent equations 




x\ = -3xi{l - xl - xl) ± \\ -yl (10) 
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xl T yl) 



yi{?>xl + 3x1 - X\ 



(11) 
(12) 



Substituting linear perturbations about the critical points xi — xi + 5xi, X2 ^ X2 + 8x2 
and Vi yi + Syi into Eqs.(fTU |) - (fT^ . to first-order in the perturbations, gives equations of 
motion 



(13) 
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5X2 


= M 


5X2 


V ^y'l ) 







where 



M 



( -3 + 9x1 + 3a;i 

6x1X2 — XVGxi 
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6x1X2 

-3 + 3xf + 9x1 - AV6x2 

6X2?/1 



txVgvi 

3x1 + 33^2 - ^Jl^i J 



\ QxiHi - Ay |?/i 

The general solution for the evolution of linear perturbations can be written as 

Sxi = Ui exp{miN) + U2 exp{m2N) + M3 exp(m3A^) 
6x2 = f 1 exp(miA^) + t>2 exp(m2A^) + fs exp(m3A^) 
6yi = Wi exp(miA^) + W2 exp(m2A^) + W3 exp^m^N) 



(14) 



(15) 



where mi, m2 and 1713 are the eigenvalues of the matrix Ai. Thus stability requires the real 
part of all eigenvalues being negative. 

S: xl + xl = l,yi = y2 = 
These kinetic-dominated solutions exist for any form of the potential, which are equivalent 
to stiff-fluid dominated evolution with a oc t^^^ irrespective of the nature of the potential, 
with the eigenvalues 



nil = 
m2 = 3 - 
ma = 6 — V6Xx2 



-Axi 
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Thus the solutions are stable to potential energy perturbations for Axi > ^/6 and Ax2 > v^. 
Using the constraint equation (|7j) we find 1 > 2xiX2 > 12/A^. That is, there exist stable 
points only for sufficiently steep (A > 2^/3) potential. 

Ai: xi 



V6' 



±(l-f),X2 = 1/2 = 
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Aa: xi = = 0, X2 = ^, 2/2 = ^±{l-^)) 
The two single-potential-kinetic solutions exist for sufficiently flat (A^ < 6) positive poten- 
tials or steep (A^ > 6) negative potentials. The power-law exponent, p = 2/A^, depends on 
the slope of the potential. From Eq. lfT^ we find the eigenvalues 

vfii = A^ 

Thus the single-potential-kinetic solutions are unstable for the positive and negative poten- 
tials. This indicates that the stability is destroyed by the potential energy perturbations of 
another scalar field. 

B: xi=x2 = yi=y2 = - S) 
The double-potential-kinetic scaling solution exist for fiat (A^ < 12) positive potentials, 
or steep (A^ > 12) negative potentials. This corresponds to a power-law solution with 
a (X t^/^\ Linear perturbations yield three eigenvalues 

mi = ^(A^-12) 

m2 = ^(A^-12)-|v/(A2-12)(A2-4/3) 
m3 = ^(A^-12) + ^V(A2-12)(A2-4/3) 

For 4/3 < A2 < 12, y/(A2 - 12)(A2 - 4/3) is replaced by zy^(A2 - 12)(4/3 - A2). So the 
terms exp(m2A^) and exp(m3A^) in the solution (jl5j) decay oscillatingly. Thus for positive 
potentials, the scaling solution is stable whenever this solution exists, whereas for negative 
potentials the scaling solution is never stable. 

The regions of (A) parameter space lead to different qualitative evolution. 

• For steep positive potentials {V > 0, A^ > 12), only a circle S exists, some kinetic- 
dominated scaling solutions of which are the late-time attractors. Thus generic solu- 
tions in the kinetic-dominated regime approach an equal-kinetic-dominated regime, 
the range of which is determined by the value of A, at late times. That is, the kinetic 
energy of each field tends to be equal via their effect on the expansion. 

• For intermediate positive potentials (V^ > 0, 6 < A^ < 12), a circle S and a fixed point 
B exist. The later is the unique late-time attractor. Thus generic solutions start in a 
kinetic-dominated regime and approach the double-kinetic-potential scaling solution. 
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Label 


Critical points 


Existence 


Eigenvalues 


Stability 


S 


xl + xl = 1, 

1/1 = 1/2 = 


all A 


(3- A^fxi); 
(6-Av%a;2);0 


stable (A^ > 12) 
unstable (A^ < 12) 


Ai, 
A2 


(^,\/±(i-¥),o,o), 


A^ < 6 (V > 0) 
>Q {V <0) 


(A^ - 6)/2; 
(A2-6)/2; A2 


unstable 


B 


1/1 = 1/2 = ^±(1 - f ) 


A^ < 12{V > 0) 
A^ > 12{V < 0) 


(A^ - 12)/4; 
|(A2 - 12)± 
|^(A2-4/3)(, 


stable {V > 0, A^ < 12) 
unstable (V < 0) 
\^ - 12) 



Table 1: The properties of the critical points 



For flat positive potentials {V > 0, A^ < 6), all critical points exist. The double- 
kinetic-potential scaling solution is the unique late-time attractor. Thus generic so- 
lutions start in a kinetic-dominated solution or in a single-kinetic-potential solution 
and approach the double-kinetic-potential scaling solution. 

For steep negative potentials {V < 0, A^ > 12), all critical points exist. The kinetic- 
dominated scaling solutions are the late-time attractors, which correspond to a con- 
tracting phase in the pre big bang scenario 



• For intermediate negative potentials (l^ < 0, 6 < A^ < 12), a circle S and two fixed 
points Ai and A2 exist. There exist no stable points. 

• For flat negative potentials {V < 0, A^ < 6), only a circle S exists, which are never 
stable. 

We now generalize the above discussion to m scalar fields, by considering each potential 
to have a different slope Aj 

Vi{(f)i) = Vojcxp {-XiKcpi) (16) 
Setting x'j = and i/^ = 0, we can get 



m 
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= 1/i ^3^x^ - Ai^Xi 

We only consider the m-kinetic-potential scaling solution where Vi ^ and Xj 7^ 0. Using 
Eq.dlHl), we get 
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Notice that XiXi/\/Q = c is an invariant. So 

m m -I 

E^| = 6c^Eti (20) 
j=i j=i '^j 

Comparing Eq.(j2ni) with Eq.jini) gives c = {'oT.T=iK'^y^ ■ From Eqs.jlTl) and (HH), we 
obtain the multi-kinetic-potential scaling solution 

- ^ (21) 



A, 




Substituting Eq. (fTIJ|) into Eq.Q gives 



(22) 



m Q 

P = Et2 (23) 



which is just the result derived in Refs.jH]. For m scalar fields with Aj = A, Eqs. fl21|) - 
yield Xi = A/(mv^), Ui = \J ±^(1 — ^) and p = 2m/A^, which are consistent with 
the above results. The multi-kinetic-potential scaling solution exists for positive potentials 
(A^ < 6m), or negative potentials (A^ > 6m). As long as each potential satisfies A^ < 2m, 
this power-law solution is infiationary. For the case m = 1, the dimensionless constant A 
must be smaller than ^/2 to guarantee power-law infiation ^2] • However, presently known 
theories yield exponential potentials with A > ^/2. In such cases multiple scalar fields may 
proceed infiation. The reason for this behavior is that while each field experiences the 
'downhill' force from its own potential, it feels the friction from all the scalar fields via their 
contribution to the expansion [S]. 

We have presented a phase-space analysis of the evolution for a spatially fiat FRW 
universe containing m scalar fields with positive or negative exponential potentials. As 
an example we study the problems of the fixed points and their stability in a two-field 
model. We find that in the expanding universe model with sufficiently fiat (A^ < 12) 
positive potentials the only power-law double-kinetic-potential scaling solution is the late- 
time attractor and the well-known infiationary solution with a (x where p = 4/A^. A 
successful infiation can be driven by multiple scalar fields with a wide range of values for 
each potential slope parameter A. We also find that the scaling solution with steep negative 
potentials is always unstable in the the expanding universe models. However, sufficiently 
steep (A^ > 12) negative potentials have kinetic-dominated solutions with a oc t^^^, which 
are always the late-time attractors. It can be known that the kinetic energy of each field 
tends to be equal via their effect on the expansion. 

We emphasize that we have assumed that there is no direct coupling between these 
exponential potentials. It is worth studying further the case their potentials have different 
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slopes. It would else be interesting to study the multi-field dynamics when a perfect fluid 
is present and for realistic cross coupling of the kind exp(Aifi;0i + A2K02)- 
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